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Abstract. - The forces acting on an optically trapped particle are usually assumed to be conserva- 
tive. However, the presence of a non-conservative component has recently been demonstrated. Here 
we propose a technique that permits one to quantify the contribution of such a non-conservative 
component. This is an extension of a standard calibration technique for optical tweezers and, 
therefore, can easily become a standard test to verify the conservative optical force assumption. 
Using this technique we have analyzed optically trapped particles of different size under different 
trapping conditions. We conclude that the non-conservative effects are effectively negligible and do 
not affect the standard calibration procedure, unless for extremely low-power trapping, far away 
from the trapping regimes usually used in experiments. 



Introduction. — The detection and measurement of forces and torques in microscopic 
systems is an important goal in many areas such as biophysics, colloidal physics and hydro- 
dynamics of small systems. Since 1993, the photonic force microscope (PFM) has become a 
standard tool to probe such forces [1-3]. A typical PFM setup comprises an optical trap - 
an highly-focused Gaussian light beam - that holds a probe - a dielectric or metallic particle 
of micrometer size - and a position sensing system. Using a PFM it has been possible to 
measure forces as small as 25 fN [4] and torques as small as 4000 fN • nm [5]. 

In order to assess the mechanical properties of microscopic systems, the first step is always 
to have an accurately calibrated optical probe. Modelling the interaction between the light 
of a focused laser beam and an extended dielectric or metallic object can be a complicated 
task [6]. The electromagnetic theory is relatively straightforward for the Rayleigh and 
the geometrical optics regimes [7,8]. However, most applications of PFM involve particles 
whose characteristic size is comparable to the wavelength of the light employed. In this 
case the exact solutions for the force-field are cumbersome to come by. Fortunately there 
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Fig. 1: (A) Schematics of the experimental setup. (B) Forces acting on a colloidal particle held by 
an optical trap. The lines represent the force-field and the color the modulus of the force. A slight 
bending of the force lines due to the non-conservative component of the force-field can be observed 
[see also Fig. 2]. 

are several straightforward methods to experimentally measure the trap parameters - i.e., 
the trap stiffness and the conversion factor between voltage and length - and, therefore, the 
force exerted by the optical tweezers on an object. The most commonly employed methods 
are the drag force method, the equipartition method, the potential analysis method and the 
power spectrum or correlation method [9]. The latter two [10,11], in particular, are usually 
considered the most reliable ones. 

An implicit assumption of all these calibration methods is that, for small displacements 
of the probe from the center of an optical trap, the restoring force is proportional to the 
displacement. Hence, an optical trap is assumed to act on the probe like a Hookeian spring 
with a fixed stiffness. This condition implicates that the the force-field produced by the 
optical forces must be conservative, excluding the possibility of a rotational component. 
This is actually true to a great extend in the plane perpendicular to the beam propagation 
direction (e.g., x-y plane in Fig. 1(B)) for a standard optical trap generated by a Gaus- 
sian beam. However, this has been shown not to be true in a plane parallel to the beam 
propagation [8,12,13] (e.g., x-z plane in Fig. 1(B)). 

Back in 1992 Ashkin already pointed out that, in principle, scattering forces in optical 
tweezers do not conserve mechanical energy, and that this could have some measurable con- 
sequences [8]. In particular, this non-conservative force would produce a dependence of the 
axial equilibrium position of a trapped micro-sphere as a function of its transverse position 
in the trapping beam (see Fig. 10(C) of Ref. [8]); such prediction was first confirmed by 
Merenda and colleagues [12]. Recently Roichman and colleagues [13] have directly investi- 
gated the non-conservative component and have discussed the implications that this might 
have for optical tweezers-based experiments making use of the thermal fiuctuations in the 
calibration procedure. 

Here we propose a technique that permits one to evaluate the relative weight of the 
non-conservative component of the optical forces. It is based on a previous work [11], where 
it was proposed an enhancement of the PFM to measure force-fields with a non-conservative 
component. We use this technique to analyze various optically trapped particles in different 
trapping conditions. The main result is that the non-conservative effects are effectively 
negligible and do not affect the standard calibration procedure, unless for extremely low- 
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power trapping, far away from the trapping regimes usually used in experiments. 

Theory. — Assuming a very low Reynolds number regime [14,15], the motion of a 
Brownian particle in the presence of an optical force-field can be described by the vectorial 
Langcvin equation 

f(i) = if(r(t)) + ViDh(i), (1) 

7 

where r{t) is the probe position and f(r) is the optical force acting on the particle, which 
depends on the position of the particle itself, of course since r is time-dependent then also 
f varies over time, 7 = Sirda is its friction coefficient, d is its diameter, a is the medium 
viscosity, ^/2D"fh{t) is a vector of independent white Gaussian random processes describing 
the Brownian forces, D = ksT/^ is the diffusion coefficient, T is the absolute temperature, 
and ks is the Boltzmann constant. 

The system that we are going to characterize is radially symmetric. It is, therefore, easier 
to study it in cylindrical coordinates. Eq.(l) can be projected in cylindrical coordinates as 



p{t) = ^fp{p,6,z) + V2Dhp{t) 
0{t) = ^fe{P^O^z) + V2D^-^ (2) 
z{t)^y,{p,e,z) + V2DK{t) 



where hp{t), hg(t) and hrlt) are independent white Gaussian random processes with unitary 
variance. 

Since an optical trap generated by a Gaussian beam is symmetrical, the particle is effec- 
tively diffusing freely with respect to the coordinate 6 and we can assume that fgi^p, 0, z) = 0. 
We can, therefore, study the movement of the particle only with respect to the coordinates 
r and z. Following a procedure similar to the one in Rcf. [11], to which we refer for more 
details, we can linearize the force-field near the equilibrium position {pq,zq) ~ (0,0) and 
rewrite the Brownian particle equations of motion as 



^p{t) + e^z{t) + V2Dhpit) 




X lu. , fKT,: ' (3) 



-e^p{t) - r^'^z{t) + ^2DKXt) 

where kp is the optical trap stiffness in the x-y plane, r]kp is the trap stiffness along the z 
axis, T] is the ratio between the trap stiffness in the x-y plane and the one along z, which is 
typically ^0.1, and e represent the relative contribution of the non-conservative component 
of the force-field. 

The terms —^p{t) and —rj^zit) represent the elastic restoring forces which are the 

conservative part of the force-field and the terms +e^z{t) and -~t^p{t) represent the non- 
conservative part of the field. In Fig. 2(A) an example of such a force-field is drawn. In 
Fig. 2(B) only the non-conservative component is depicted. 

The statistical analysis of the Brownian trajectories permits one to reconstruct the force- 
field acting on the particle. Here we will use the autocorrelation functions (ACFs) of the 
radial (p) and axial (z) particle position, and the difference between the two cross-correlation 
functions (DCCF) between the the radial (p) and axial (z) particle position. 

Assuming e << 1 and e << 77, the ACFs for p and z are decoupled 

ACF,(r) = :^exp(^-^|r|) (4) 

and 

7^ ( Vkp I 



ACF,(r) = ^cxp(^-^|r|). (5) 
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Fig. 2: (A) Force field generated by an optical trap in the presence of a rotational component and 
(B) the rotational part of the force-field with rj = 0.3 and e = 0.1. Note that these values are 
larger than typical ones, which are usually about rj ~ 0.1 and e ~ 0.05, in order to have a clearer 
qualitative picture of the way they affect the force-field. 



These expressions are the standard ones for a conservative force-field and are independent 
from e [11]. The DCCF is 



(6) 



sinh 



v/|(l-ry)2-4e2| 



Equations (4), (5) and (6) can be used to fit the values for fcp, ij, and e, and therefore to 
reconstruct the force-field up to the first order. In particular, the ACFs (4) and (5) can be 
used directly to fit the values of kp and rj. Once these values are estimated, it is possible to 
fit the value e by using the slope of the DCCF (6) around r = 0. 

Once the force-field parameters have been fitted, it is possible to calculate the torque 
acting on the particle due to the presence of the rotational component of the force-field [5,11] 
as 

T = ekp [Var(p) + Var(z)] , (7) 

and the circulation rate as 

n = e-^. (8) 
27r7 ^ ' 

Experimental setup. — The experimental setup, shown in Fig. 1(A), is described 
in details in Ref. [16]. The PFM comprises a home-made optical microscope with a high- 
numcrical-apcrturc water-immersion objective lens (Olympus, UPLAPO60XW3, NA=1.2) 
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and a frequency and amplitude stabilized Nd-YAG laser (A = 1.064 /im, 500 mW maximum 
output power, Innolight Mcphisto). 

Polystyrene micro-spheres (Scrva Electrophoresis, 1.06g/cm'^ density, 1.65 refractive in- 
dex) with a diameter of 0.45 ± 0.01 /xm and 1.25 ± 0.05 /xm were diluted in distilled wa- 
ter to a final concentration of a few particles//il. A droplet (100^1) of such solution was 
placed between a 150 /xm-thick coverslip and a microscope slide, which were separated by 
a 100 ^m-thick parafilm spacer and sealed with vacuum grease to prevent evaporation and 
contamination. Such sample cell was mounted on a closed-loop piezoelectric stage (Physik 
Instrumente PI-517.3CL), which allowed movements with nanometer resolution. The sam- 
ple temperature was continuously monitored using a calibrated NTC thermistor positioned 
on the top surface of the microscope slide and remained constant within 0.2 degrees during 
each complete set of measurements. 

A micro-sphere was trapped and positioned in the middle of the sample cell, i.e. far 
away from the glass surfaces to avoid hydrodynamic effects on the bead motion [10]. Its 
3D position was monitored through the forward scattered light imaged on a a InGaAs 
Quadrant Photodiode (QPD, Hamamatsu G6849) at the back focal plane of the condenser 
lens [17], using a digital oscilloscope (Tektronix TDS5034B) for data-acquisition. The QPD- 
response was linear for displacements up to 300 nm (2nm resolution, 250kHz bandwidth). 
The conversion factor from voltage to distance was calibrated using the power spectral 
density method [18]. We excluded significant deviations from a harmonic trapping profile 
by verifying that the trapping potential could be fitted to such a profile in all the three 
directions. 

Experimental results. — We performed the experiments using particles with diameter 
0.45 ^m and 1.25 //m. Similar particles are commonly employed in experiments that use 
optical traps [3]. The particle positions were acquired at 2.5 kHz, which is above the cutoff 
frequency of the particle motion in the optical trap and below the QPD bandwidth. 

For a dataset of 2N-I-1 particle position - i.e. a;„, y„ and 2„ for n = —N, —1,0, 1, TV 
at times tn = 0.4 • n ms - the experimental ACFs and DCCF are 

N 

ACF(,^)(T = 0.4-m) = J2 P^n-nPn, (9) 

N 

ACF(^)(T = 0.4-m) = ^ z^.nZn. (10) 

and 

N 

DCCF(,^)(r = 0.4-m) = ^ p„_„z„ - z„_„p„, (11) 

n=-N 

where p„ = y^x"^ + y'^ and z„ are the time-series of the particle position in cylindrical 
coordinates. 

For medium-high optical power at the sample (few tens of milliwatts) the DCCF|,'^^ (t) = 
within the experimental error. This is an evidence that the contribution of the non- 
conservative force-field component is effectively negligible, or at least undetectable for ac- 
quisition times up to several tens of minutes. 

In order to have a non-vanishing DCCF the optical power at the sample needed to 
be reduced down to a few milliwatts. Furthermore, even under such low power a clearly 
non-vanishing DCCF was only obtained when the acquisition time was increased up to 
400 s (i.e., 2Af 1 w 1 X 10^ particle positions). It is an important remark the fact that, 
due to extremely low value of the rotational contribution to the total force-field, it is not 
immediately evident from the traces of the particle motion. Indeed the particle undergoes a 
random movement in the p — z plane, where it is not possible to distinguish the presence of 
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Fig. 3: Autocorrelation functions (ACFs) and difference of tlie cross-correlation function (DCCF) 
for a 0.45 ^.m colloidal particle optically trapped with a laser power of 6 mW at sample (first line of 
Tab. 1). Inset: close-up of the DCCF. Dashed lines represent experimental data, while solid lines 
are the curves obtained from the fit with Eq.(4), (5) and (6) 

a rotational component without the aid of a statistical analysis such as the one we propose 
(see videos in the supplementary materials). 

In Fig. 3 the ACFs and DCCF are presented for a 0.45 ^m diameter particle held in an 
optical trap with an optical power at the sample of 6.0 mW (average of 10 400 s-datasets). 
A good agreement is found between the experimental and theoretical ACFs and DCCF. In 
particular, the experimental and theoretical DCCF are very similar over all the range of 
T even though the fitting was performed only on the central slope. The amplitude of the 
DCCF is very small compared to the ACF one and the stiffness of the optical trap is quite 
low, only 5 pN/ fim, compared to common experiments with optical tweezers. 

In a very weak optical trap the particle can explore regions far away the trap-center. 
Thus the particle motion is more influenced by the non-conservative force. This means that 
the less power is used the more the DCCF amplitude is observed. This is shown in Figs. 4(A) 
where we can see the behavior of the DCCF for a 0.45 /im diameter particle while decreasing 
the optical power. The DCCF amplitude increases as the power decreases. Furthermore, 
the range over which the DCCF is not vanishing broadens. This translates into an increase 
of the rotational component relative weight for decreasing power as it is shown in Fig. 5(B). 

In Fig. 4(B) we can see similar traces for a 1.25 fim. Again there is a good agreement 
between the experimental and theoretical DCCF over a wide range, even though the fitting 
was performed only on the central slope. Again the rotational component relative weight 
increases for decreasing power (Fig. 5(B)). 

It is worth to be noted that if we compare the DCCF curves for the two diameters used 
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Fig. 4; Cross-correlation function for an optically trapped particle whose diameter is (A) 0.45 fim 
and (B) 1.25 fim. The powers used are listed in Tab. 1 and the color legend is: P(blue) > P(red) > 
P(green). Dashed lines represent experimental data, while solid lines are the curves obtaining from 
the fit 



at comparable stiffness, i.e. kp = 2.9 pN/fim for d = 0.45 /im and kp = 3.6 pN/nm for 
d = 1.25 /im (second and last line of Tab. 1) we can observe a slightly larger amplitude for 
the smaller particle. Again this is a confirmation that the effect is due to the larger volume 
explored by the smaller particle with respect to the larger one. 

In the Tab. 1 the numerical values for the studied cases are presented. The stiffness 
along the horizontal plane kp is also shown in Fig. 5(A) where it can be appreciated the 
fact that it is linear with respect to the optical power. The ratio of the stiffness along z and 
p does not depend on the power used, while a clear increase of the rotational component 
e is observed in both particle diameters used in this work. We notice that such values are 
similar to the ones reported in Ref. [13]. In their case the trap power is larger than the one 
used in the present experiment, but for the different setup the resulting stiffness is lower 
due to the fact that the particle size is larger (2 //to). Nevertheless the order of magnitude 
of the circulation they observed is in agreement with our data. 

In particular the torque associated to the non-conservative component of the force-field 
is calculated according to Equation 7. The resulting values are extremely small (see Tab. 
1), being actually orders of magnitude smaller than the ones previously reported, e.g., 4 x 
10^ fN ■ nm for the torque transfer from a Laguerre-Gaussian beam to a Brownian particle 
[5], 6 X 10^/7V • nm for microscopic hydrodynamic flows [19], 1 x lO^/A^ • nm for DNA 
twist elasticity [20], 5 x 10^ fN ■ nm for the movement of bacterial flagellar motors [21], 
2 X 10^/7V • nm for the transfer of orbital optical angular momentum [22] , or 5 x 10^/7V • nm 
for the transfer of spin optical angular momentum [23]. However, we must remark that in 
this letter we are focusing on the non-conservative forces that arise in a standard optical trap 
due to the fact that the trap is not perfectly harmonic in the vertical plane, while some of 
the previously mentioned examples, such as experiments with light beams that carry orbital 
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Table 1; The experimental parameters and the values obtained from the fit are reported here. 
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z— stiffness 
p— stiffness 


rotational 
component 


torque 


circulation 
rate 


d inm) 


P (mW) 


kp (pN/^m) 


V 


e 


T {fN ■ nm) 


n (Hz) 


0.45±0.01 


6.0±0.2 


5.0±0.1 
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l.OiO.l % 
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2.1±0.3 
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1.7±0.2 % 
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1.68±0.07 
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2.4±0.2 % 
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1.25±0.05 


4.0±0.2 


14.0±0.7 


14±1% 


0.7±0.1 % 


210±30 


1.4±0.2 


1.25±0.05 


1.4±0.2 


4.4±0.2 


13±1% 


1.2±0.1 % 


420±70 


0.8±0.1 


1.25±0.05 


1.0±0.2 


3.6±0.2 


12±1% 


3.8±0.3 % 


1400±160 


2.1±0.3 



or spin angular momentum, refer to rather different situations in which the dominant effect 
is the presence of a non-conservative (or rotational) force field, which generates the particle 
movement. 

Conclusions. — The DCCF^'^],(r) = within the experimental error is a clear evidence 
that the contribution of the non-conservative force-field component is effectively negligible. 
This is actually true in most practical experimental situations. Typical optical tweezers 
experiments are performed with various tens of milliwatts of optical power at the sample 
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Fig. 5: (A) Behavior of the radial stiffness kp as function of the laser power for the two kinds 
of particle used. The solid lines are the linear fits. (B) Value of the relative contribution of the 
non-conservative component e as a function of the power used. Notice that it increases for lower 
powers. 
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and acquiring data for at most a few minutes [3]; however, we needed to decrease the optical 
power at the sample down to a range of a few milliwatts and to increase the acquisition 
time up to several minutes in order to see the signature of a non-conservative force-field 
component in the DCCF. 

In particular, for medium-high laser powers (from a few tens of milliwatts at the sample 
onwards) the effect is undetectable even for long acquisition times - indeed, the value of e 
steadily decreases as the laser power is increased (see Tab. 1 and Fig. 5(B)). This is due to 
the fact that the deviation from a conservative force-field is larger far from the trap center, 
which is explored more often in a weaker trap. Furthermore, this can also be seen in Figs. 
6 and 7 of Ref. [12]. 

Therefore, we conclude that for the trapping regimes that are usually employed in ex- 
periments, the effect of the non-conservative force-field component are effectively negligible 
and do not affect the standard calibration procedure. Whenever a doubt is present, the 
extension to the standard optical tweezers calibration techniques we have proposed in this 
letter can be used to verify to what extend the assumption of a conservative force-field is 
fulfilled. In particular this technique might become increasingly useful as optical measure- 
ments of forces reach to ever smaller length and force scales, because these nonequilibrium 
effects will become increasingly noticeable. 

* ^ * 
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